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Abstract. We study the edge deletion process of random graphs near a fc-core 
percolation point. We find that the time-dependent number of edges in the process 
exhibits critically divergent fiuctuations. We first show theoretically that the fc-core 
^ ' percolation point is exactly given as the saddle-node bifurcation point in a dynamical 

^ ■ system. We then determine all the exponents for the divergence based on a universal 

description of fiuctuations near the saddle-node bifurcation. 

> 

^ ; PACS numbers: 05.10.Gg,05.70,Jk, 64.60.ah 

'-^ . 1. Introduction 

Q I We study the following time evolution of random graphs on n vertices. Let us denote one 
sample trajectory of graphs by {G{t)}o<t<oo, where the time t is a real number. G(0) is 
I assumed to contain m edges that connect two vertices chosen randomly. Here, R = m/n 
\^ I is regarded as a control parameter of the model. When t > 0, a vertex is chosen with a 
^ ' constant rate a for each vertex. Then, if the number is less than an integer k, all the 
O I edges incident to the vertex are deleted. This rule defines the Poisson jump process G(t) 
00 ■ in the set of graphs on n vertices. We display two examples of the time evolution of 
(yQ ' graphs in Fig. [H where the random graphs are embedded in the two-dimensional space. 
I Let be the number of edges at time t. Obviously, when R is sufficiently small, 

fi{t = oo)/n is zero with probability 1 in the limit n — >■ cx). It has been known that 
^ ' there is a critical value of R above which /x(oo)/n is finite (nonzero) with probability 1 
■ in the same limit [I]. The final graph provides the fc-core (See Fig. [1]), which is defined 
by the largest subgraph with minimum degree at least k. This transition with respect 
to the change in R is called the /c-core percolation in a random graph. The critical value 
Rc was calculated exactly in Ref. [I]. 

The fc-core percolation was studied in several research fields such as magnetism [2], 
rigidity percolation [3], jamming transitions [U [5], and network problems [HI IZ]- It is 
also related to the random field Ising model, which is a representative model exhibiting 
so-called avalanches [HI E] . In particular, the dynamics of fc-core percolation might be 
considered from the viewpoint of the vulnerability of a network to random node attack 

§ To whom correspondence should be addressed (sasa@jiro.c.u-tokyo. ac.jp). 
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Figure 1. Time evolution of random graphs from initial states (left) to final states 
(right), k = 3, n — 118, and R = 1.0 (a) and R = 1.7 (b), respectively. The final state 
in the case R = 1.7 corresponds to the 3-core. The filled square symbol (red online) 
at each time represents a heavy vertex and open square symbol (blue online) a light 
vertex. 

In this paper, we wish to elucidate the nature of the dynamics near the transition 
point. Concretely, let h be the number of vertices with degree at least k. (Such a vertex 
is called a heavy vertex; otherwise, a light vertex.) We are interested in the time evolution 
of h. As an example, we present the results of numerical simulations in Fig. [2l [Jj] Here, 
the ensemble average {h{t)/n) and its fluctuation intensity x(t) = {{h{t) — {h{t))Y/n) 
are displayed as functions of t. Figure [2] indicates that x(t) has one peak at t = r, and we 
conjecture that r and xi'^) exhibit the power-law divergences r ^ and x{t) — ^ 
where e = -Rc ~ -R > 0. Indeed, we will derive these divergences theoretically and 
determine the values C = 1/2 and 7 = 5/2. 

The divergent behavior observed near the percolation point suggests the existence 
of critical fluctuations. On the other hand, it has been known that a giant fc-core appears 
in the discontinuous manner at the transition point for cases /c > 3. Such coexistence 
of the discontinuous transition and critical fluctuations has been emphasized in relation 
to the nature of jamming and glass transitions [HI [121 US]- Therefore, the theoretical 
description of the divergent behavior near the fc-core percolation may provide a new 
insight toward understanding of jamming and glassy systems. 

This paper is organized as follows. In Sec. [21 we present a master equation for four 
variables that characterize a graph. Since this master equation was derived rigorously in 
Ref . [1] , our presentation in this paper is based on an intuitive argument understandable 
for physicists. Then, in Sec. [3l by considering the situation with large n, we derive a 
Langevin equation for the four variables. The Langevin equation is analyzed in the 
subsequent two sections. In Sec. [H we find a saddle-node bifurcation for the rate 

II In numerical simulations, we generate a chain of waiting time obeying the Poisson distribution and 
choose a vertex randomly with these time intervals. 
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Figure 2. Relaxation behavior of heavy vertex density {h{t)/n) (a) and its fluctuation 
intensity xW (b). e = 0.03, e = 0.05, and e = 0.07. /c = 3, a = 1, and n = 4096. 

equation obtained by the limit n oo. Here, the bifurcation point corresponds to the 
/c-core percolation point. Then, in Sec. [5l we study effects of noise near the bifurcation 
point and calculate the exponents that characterize critical divergences. The final section 
is devoted to concluding remarks. In order to simplify the argument, we consider the 
case k = 3. The generalization to cases > 3 is straightforward, and essentially the 
same results are obtained. 

2. Master equation 

Let At be a sufficiently small time interval. We can describe the stochastic process by 
the transition probability P{G'\G), which is the probability that G{t + At) = G' under 
the condition that G{t) = G. Since P{G'\G) is a huge matrix, we cannot treat it directly. 
Hence, we wish to have a simple description of the dynamics. The simplification of the 
dynamics of G{t) consists of two steps. In the first step, we describe the dynamics in 
terms of the characteristic quantities of the graph such as the number of edges fj, and 
the number of vertices Vr with degree r, where r = 0, 1, ■ ■ ■. Among them, the number of 
light vertices, vq, vi and V2 are directly related to the dynamics of the graph because all 
the edges incident to a chosen light vertex will be deleted in the next change of the graph. 
Indeed, according to Ref. [1], \^ the time evolution of the four-tuple w = (/i, fo, "Wi, ^2) 
is described by a Markov process. Mathematically, the probability of w' at time t + At 
provided that w at time t is given, which is denoted by p{w'\w), is expressed as a 
function of w' and w for general n. (See proposition 1 in Ref. [1].) Subsequently, in 
the second step of the simplification, the asymptotic formula of p{w'\w) for large n is 
derived. (See Corollary 1 in Ref. p.) 

In this paper, we do not review the derivation of the asymptotic form p{w'\w) in 
Ref. [1]. Instead, we provide its mathematically naive derivation by focusing on cases 

% Note that the edge deletion process in Ref. [P is not identical to the dynamics we define. First, their 
dynamics are given as a discrete process. Second, in their dynamics, a non-isolated light vertex is always 
chosen at each time step. Thus, as time goes on, the deletion in their dynamics is accelerated more 
than that in our evolution rule. Despite this difference, one can transform mathematical statements in 
[1] to those valid in our model. 
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with large n from the outset. More precisely, we estimate p{w'\w) assuming that Vr/n 
with r > 3 takes the most probable value hq^/n, where from the law of large numbers, 
qr is equal to the probability that r edges are incident to a given heavy vertex under the 
condition that w is specified. (Recall that h represents the number of heavy vertices 
that is equal to J2r>k=3'^r- ) We express this statement formally as 

J = 

where the probability qr is given by the Poisson distribution 

with the normalization constant 

Q^z) =e-'{e' -1- z- zy2). (3) 

Although the appearance of the Poisson distribution seems natural, its mathematical 
proof is not simple. (See Ref. [1].) Here, using the trivial relatiorCJ 2/i — t>i — 2t>2 = 
^^grfr, we have 



h 



5^rg„ (4) 

r=3 

which determines z in ([2]) and ([3]) for a given w. Since the direct calculation using ([2]) 
leads to 

with 

U{z) =e-'{e' -1- z), (6) 

we obtain a useful relation for the determination of z from w: 

2f 2 _ zli{z) 

h ~ Q{z) ■ ^ ' 

In the argument below, z always represents the unique solution of ([7]) for a given w. 
Furthermore, one can easily confirm the relation 

z = ^rqr, (8) 

which provides us a simple interpretation of z. 

Now, we estimate p{w'\w). We first notice the value of /i' — yU. (Note that /i' — 
represents the change of the number of edges during the time interval At.) (i) When 
/i' — /i = 0, no deletion occurs. This implies w' = w. (ii) When /i' — /i = —1, one 
edge incident to a chosen vertex is deleted. This edge connects the chosen vertex with 

+ We express an edge by a pair of integers (^1,^2) when the edge hnks a vertex of degree £1 and 
another vertex of degree £2- Collecting all the pairs of integers for the edges, we have 2/x integers. 
Here, count the number of integers that are greater than 2. 
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(-1 
(-1 
(-1 
(-1 



,2,-2,0) 
,1,0,-1) 
,2,-1,1) 
,1,-1,0) 



avipi/n 
aviP2/n 
avips/n 
aviPi/n 



Table 1. Rate r{(Tj\w) for the transition w ^ w + (Tj. pj = jvj/{2ii) for 1 < j < 3, 
and p4 = zh/{2ii). 



another vertex with degree £. Then, w' — w takes four values depending on £ = 1, 2, 3, 
and i > 4:, which are denoted by crj with j = 1, 2, 3 and 0-4, respectively. (See Table. 
[U) (iii) When fi' — fi = —2, two edges incident to a chosen vertex are deleted. Each edge 
connects the chosen vertex to another vertex with degree it, i = 1,2, where ii = 1, 2, 3 
or ii > 4. We assume ii < I2 without loss of generality. Then, w' — w takes ten values 
depending on the values of ii and £2, which are denoted by o-j, 5 < j < 14, where 
the correspondence between j and (£1,^2) is shown in Table O (iv) We do not need 
to consider the cases /x' — yU < —3. Although such cases appear when deletions occur 
twice or more during the time interval, the probability of their occurrence is negligible 
for sufficiently small At. To sum up, w' — w takes either or crj, j = 1, ■ ■ ■ , 14, and 
the fourteen transitions occur independently. 

We denote the rate of transition w ^ w + cr hj r{cr\w). We can then write 



when w' ^ w. 6{x,y) is the four- dimensional Kronecker delta function ior x, y in N^. 
p{w\w) is determined from the normalization condition of the probability. 

Let us estimate the transition rate r{crj\w). Recall that a vertex is randomly chosen 
at the rate a (per unit time and per each vertex). Then, the probability that the degree 
of the chosen vertex equals to one is given by vi/n. We next consider the conditional 
probability that the edge incident to the chosen vertex connects it with a vertex of 
degree r. This probability, which is denoted by c^, takes a complicated form for general 
cases. (See Ref. [T]). Here, notice that rVr/^n is the probability of finding a vertex with 
degree r when we observe one vertex connected to an edge, which we choose randomly. 
The difference between and rfr/2/i originates from the condition under which an edge 
is chosen. The difference is negligible for sufficiency large n. Combining these results, 
the rate r{aj\w) is estimated as 



14 




(9) 



r{(Tj\w) 




2/i 



(10) 



for 1 < j < 3, and 



00 




(11) 
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10 


(2,3) 


(-2,1,1,-1) 


2av2P2P3/n 


11 


(2, > 4) 


(-2,1,1,-2) 


2av2P2PA/n 


12 


(3,3) 


(-2,1,0,1) 


av2p1/n 


13 


(3, > 4) 


(-2,1,0,0) 


2av2P3P4/n 


14 


(>4,>4) 


(-2,1,0,-1) 


av2pl/n 



Table 2. Rate r{(Tj\w) for the transition w ^ w + a-j. pj = jvj/{2^) for 1 < j < 3, 
and p4 — zh/(2fi). 



Here, it should be noted that V3 and ^^4^ rvr in flTUl) and ffTTl) are calculated from 
([1]) with z determined by ([7]). For convenience of later calculation, we summarize the 
result in Table [U In this table, we introduce pj = jvj/{2fi) for 1 < j < 3, and 
Pi = ?''f r/(2yu) = zh/{2fi). (We used ([8]) in deriving the latter equality.) The 

transition rate r{(Tj\w) with 5 < j < 14 is calculated in the same manner by noting 
that two edges incident to one vertex can be treated independently. The result is 
summarized in Table [2l In the argument below, we set a = 1 for simplicity. 

Before closing this section, we consider the initial condition w{0). In order to 
simplify the argument, we assume that w{0)/n takes the most probable value in the 
limit n — > 00. Let us calculate this value. We first consider the probability that k edges 
are incident to a vertex chosen randomly: 

p {n-l)Ck ■ {N-n+l)C(rn-k) „s 

- , l-L^j 

where = „C2. Taking the limit n 00 with fixing R{= m/n), we obtain 

= ^(2i?)'=e-2^. (13) 
kl 

This leads to 

w{0) = (m, n{2R)e'^'^, n{2Rfe-^''/2). (14) 

3. Langevin equation 

We define the density variable p = (/i, po. Pi, P2) by w/n. When n is large, the dynamics 
of p are expected to be described by a Langevin equation. We shall derive the equation 
from the transition probability p{w'\'w) given in ([9]) with Tables [1] and [2l 
We utilize an expansion formula of the Kronecker delta function 

5(^,2/)= r^e-(--^) (15) 
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for any x and y in N. Substituting this formula into ([2]), we obtain 

p{np'\np) =nAt -— ^e^'^ ^'''-'') ^ r(cr, |np)e'"-"^/" (16) 



^14 



for p' 7^ p, where 2? = [—n7i,mi]^. Noting p{np\np) = 1 — r(crj |i(;)nAt, we can 

write 

for any p' and p, where JF is expressed as 

14 / 14 

J^(p) = riAt ^r(o-j|np)e^^-'^^/" + I 1 - ^ r(<Tj|np)nAt 

i=i V i=i 

/ 14 

= exp J]r(o-j|np)nAt (e*" '^^/" - l) + 0{/\t^ 
\j=i 

/ 14 \ 

~ exp ^ r(crj|np)nAt (in • CTj/n — (tt ■ crj)^/(2n^)) j , (18) 

where we have ignored the terms of 0((At)^, At/n^). Here, we first consider the case 
with large n and then assume At to be sufficiently small. We then have the transition 
probability 

/ 14 

—iu ■ j p' — p — r{(Tj\np)At(Tj 



p{np'\np) = I — — rjexp 



R4 " 



(27m)4 



14 ^ ^ 



(19) 



Here, we define a 4 x 4 matrix 

^.n.(np) = '-^^^{<tM<^,)^. (20) 

Since the matrix A = (Aim) is semi-positive, there exists the semi-positive matrix G 
satisfying A = G^. We also define AH by 

14 

G(np)AH = p' - p - ^r{(Tj\np)At(Tj. (21) 

Then, the probability density of AH is expressed as 

p(AH) = det(G)nV(np'|np), (22) 

where det(G) is the determinant of the Jacobian matrix associated with the 
transformation from p' to AH, and note that the probability density of p' is given 
by n^p{np'\np), because j d'^ p'p{np'\np) = 'Y^^, p{w'\np) / n'^ = From (|T9|) and 

221), we obtain 
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This implies that AS is the Gaussian noise satisfying 

2Af 

(AS^AS^) = 6i^. (24) 

n 

Taking the hmit At -> in ([21]) with we obtain 

dtp = '^r{(Tj\np)(Tj + J-G{np) ■ ^ (25) 
i=i 

with {C,i(t)C,m(t')) = 26(t — t')5im- Here, the symbol " ■ " in ( l25l) represents the Ito rule 
of the multiplication of stochastic variables. Finally, from (fT^ . the initial condition of 
the Langevin equation is given by 

p(0) = (i?, e-2«, (2i?)e-2^, (2i?)2e-2^/2). (26) 
4. Deterministic equation 

The Langevin equation fl25l) becomes the deterministic equation in the limit n ^ oo: 

14 

<9tp= ^r(crj|np)o-^-. (27) 

Concretely, using Tables [1] and [21 we can obtain the expression 

dtji = -s, (28) 

dtPo = Pi + P2 + Pi—, (29) 
2/i 

dtPi= -Pi-Pi^ + 2P2^, (30) 

9,p2= -P2-2p2^ + 3p3|r, (31) 

where s = pi + 2p2 is the density of the edges incident to light vertices, and ps = v^/n 
is determined as a function of p from ([1]) with z determined by ([7[). The derivation 
of (l28 l) -( l3Ti) requires tedious calculation, while the result is understood intuitively. For 
example, the third term of ([29l) represents the change in the degree of a vertex from 1 to 
by the deletion of the edge connecting this vertex with another vertex that is chosen 
randomly. The differential equation given in f[28|) - f[3T|) with the initial condition f[26l) 
determines the most probable behavior of the Langevin equation ( l25l) . 

Putting aside the initial condition, we study the differential equation in ( l28i) - ( 1311) . 
The complicated nature arises from the implicit dependence of pa on p through z. In 
order to avoid it, we carry out the transformation of variables. First, we choose z as a 
dynamical variable. Taking the derivative of ([7]) with respect to time, we obtain 

(32) 

Seeing f[28|l and (!32|) . we further choose s as a dynamical variable. We then obtain 
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Thus, (|5UI) . ([32D, and constitute the differential equation for {fi, pi, s, z), 

which is equivalent to the differential equation for p. Here, from fl28l) and fl32|) . we 
find immediately a constant of motion 

^2 



(34) 



Furthermore, noting dth = —3p3s/{2fi) and dzQ = z^e ^/2, one confirms that there is 
another constant of motion 

h 



J,. 



(35) 



where h = h/n. Recalling ([7]), we rewrite J2 as J2 = zll{z)/{2p — s). Now, defining 
-^3 = J1J2/2, we obtain the expression 

U(z)\ ^ 
Ji 



s = 2\l- J- 



(36) 



which defines integral curves in the (2;, s) space. 

Using the constants of motion, we choose a set of dynamical variables as ^ = 
{Ji, J3, Pi, z). Then, the differential equation for ( Ji, J3, 2;) takes the simplest expression 
that dtJi = 0, dtJs = 0, and 

dtZ = -z + JsU{z). (37) 

Importantly, the time evolution of z is independent of pi. Besides, the initial condition 
(!26l) leads to Ji = 4i? and J2 = 1; thus, J3 = 2R. To sum up, the dynamical behavior 
of the edge deletion process of random graphs is described by 

dtz = -z + 2RU{z) (38) 

with z{0) = 2R. 





Figure 3. (a) The graphs of z and 2i?n(z) with R = 1.675. Two intersection points 
appear when R is increased, while no intersection when R is decreased, (b) Shapes of 
potential functions V{z) for the cases that R ~ Rc, R < Rc and R > Rc, respectively. 

The differential equation given in ( !38l) can be easily analyzed. First, there exists 
the trivial solution z = 0. Then, let z^, be another fixed point (if it exists). -2^(7^ 0) 
satisfies z^ = 2RIl{z^,). From Fig. [3l we find that two nontrivial solutions exist when 



Dynamics of k- core percolation in a random graph 



10 



R > Rc, where Rc is determined by min2>o[2^ ~ '^RJ^iz)] = 0, which yields the compact 
expression of Rc as 

i?c = -min ( ) . (39) 

We numerically calculated -Rc — 1-675. In order to investigate the solution trajectory of 
(1371) . we express (1371) in the form 

dV 

4. = -^. (40) 

where the potential function V{z) is given by V{z) = z'^/2 — 2R{z + 2e~^ + ze~^). 
We display the shape of the potential V{z) in Fig. [31 It is seen that there are a 
pair of minimum and maximum in addition to the trivial minimum point z = when 
R > Rc- Obviously, the solution corresponding to the maximum (saddle) is unstable, 
while the solution corresponding to the minimum (node) is stable. Note that the 
potential is a monotonic increasing function in z when R < R^, which corresponds 
to the fact that there is no nontrivial stationary solution when R < R^. The qualitative 
change of trajectories ai R = Rc is called a saddle-node bifurcation. The fixed point at 
R = Rc is the marginal saddle, which is denoted by Zc{= z^,{Rc)). Since the condition 
z(t oo) 7^ represents the existence of a fc-core, Rc is the fc-core percolation point. 
This determination method of Rc is essentially equivalent to that in Ref. |T]. The 
achievement of this study is the identification of the bifurcation type observed in the 
dynamics of fc-core percolation in a random graph. 

Now, we investigate the behavior of the system with R = Rc — e, where e is a small 
positive constant. We define a dynamical variable (p hj z = Zc + (p- Substituting it into 
( !38ll . we obtain 

dt(j)= -ea-b(j)^ + 0{(j)^), (41) 

where a = 2Il{zc) and b = —Rcll"{zc). Since the solution 0(t) is expressed as a scaling 
form 

(f){t) = ei/'$(e^/'t), (42) 
the typical time for exiting the marginal saddle is proportional to e"^/^. 



5. Critical fluctuation 

Next, we study the fluctuations that are described by the Langevin equation (!25|) . Since 
the deterministic equation for ^ is the simplest one, we rewrite the Langevin equation 
by using ^. Formally, we express the variable transformation from p to ^ as ^ = (pip)- 
Then, using Ito's formula, one can derive 

a,. = + J,n(.) + X: ^^-G,,?, + - Y: (43) 

where the expression dz/dpj is evaluated from the functional dependence of z on p. It 
should be noted that J3 fluctuates in the Langevin description. 
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Concretely, we investigate the divergent behavior of the quantity 

X.{t) = n{{z{tr) - {z{t)f) (44) 

for the system with R = R^ — e, where e is a small positive constant. Here, it is naturally 
expected that the divergent part of the fluctuations of h{t) is identical to that of z{t). 
We therefore conjecture that Xz{t) has a peak at t = and that and Xzit^,) exhibit 
the power-law divergences ^ e^^ and Xz{t*) — ■ We shall derive these divergences 
theoretically. 

The perturbative calculation with respect to the nonlinearity in fH3|) seems quite 
difficult to capture the divergent behavior of z{t). Instead, we utilize the bifurcation 
structure, as done in Ref. [ll[J. Following the idea of the method, we first notice two 
solutions of (138!) with R = R^. One solution satisfies the conditions z-^it) Zc for 
t — oo, Zuit) — > for t ^ oo, and z^^{0) = zq, (say, zq = 0.5). The other solution 
ZB{t) satisfies the conditions ZBit) ^ Zc ior t ^ oo and -2b(0) = -2(0). zb(^) represents 
the down-hill trajectory to the marginal saddle from the point Zo in the potential shape. 
Then, we express the trajectories by using the exit time 6 from the marginal saddle as 

z{t) = z^it ~e) + {zB{t) - Zc) + ^{t - e), (45) 

where (p{t — 6) represents a deviation from the superposition of the two solutions. 

It is worthwhile to note that the variable 6 corresponds to the Goldstone mode 
associated with the time-translational symmetry. Thus, the fiuctuation of 6 carries 
a divergent part, while can be treated as a variable slaved to 6. Based on this 
observation, {z{t)) and Xz{t) can be estimated by the statistical average over 6. Note 
that we have devised a theoretical framework in which the statistical distribution of 6 
can be calculated perturbatively by considering the interaction of 9 with {zB{t) — 1) 
[m Us]. In the argument below, without entering this lengthy calculation, we shall 
determine phenomenologically the exponents characterizing the divergent behavior of 
Xz{t). 

We first calculate the exponents characterizing the divergences of {6) and the 
intensity of fiuctuation xe defined by 

^,^ni{9')-{9f). (46) 

We start with the scaling relations 

{6) =n^'/-/,(ni/-e), (47) 

xe = n^'/-/2(ni/'^-e), (48) 

for large n and small e, where we have introduced the exponents v^, and 7'. We here 
assume that /i(0) = const and /2(0) = const. We also assume that /i(x) ~ x~'^' and 
/2(x) ~ x~'~^' for X ^ 1, because {0) and xe are expected to be independent of n in the 
regime x ^ 1. We find that C,' = 1/2 from (H2ll . Furthermore, from (H7j) . we assume 
that a distribution function of 9 is expressed as an ra-independent function of 9n~'^^'^* 
when e = 0. This leads to a relation 'j'/i'* = 2(^'/i/^, + 1, which yields 

7' = 2C' + z/*. (49) 
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To this point, we have avoided the analysis of (1431) . In order to determine the value 
of u^, we need to study the equation. For the sake of a simple argument, we assume that 
the behavior near the saddle-node bifurcation point is described by ( HTl) with a noise 
term: 

dt(l) = -ea - + \ -ri, (50) 
V n 

where rj satisfies {ri{t)ri{t')) = 6{t — t'), and is a constant. Here, we have ignored 
effects of fluctuations of J3 and the variable dependence of noise intensity. It is a non- 
trivial mathematical problem to clarify whether these simplifications are allowed in the 
description of critically divergent fluctuations. Q 

Once we are allowed to use (|50l) . we can determine the value of the exponent z/* as 
follows. We set e = and write the weight for trajectories [0] = (0(t))o<t<oo: 



P([0]) = |exp 



2a n 



(51) 



where the last term originates from the Jacobian term associated with the 
transformation from rj to (p. Now, we define a new scaled variable $(s) by (j){t) = 
n~^^^^{s) with a scaled time s = n~^^^t. Substituting this into ( ISTl) . we can confirm 
that the distribution function of trajectories $(s) is independent of n. This implies that 
the time scale near the marginal saddle is proportional to n^^^. This yields C/^* = 1/3. 
Recalling (' = 1/2, we have arrived at z/* = 3/2. From (H9ll . we also obtain 7' = 5/2. 
The result is summarized as follows. 

(6) =n'r'f\{n'/\), (52) 
Xe =n'/'Mn'/h). (53) 

From these, (9) ~ e^^/^ and xe — e^^/^ in the regime 0(n~2/^) ^ e ^ 0{n^). It should 
be noted that (|52|) and fl53|) have been confirmed numerically for a simple stochastic 
differential equation whose local form near the marginal saddle is equivalent to ( !50l) [15j . 

Now, using this result, we calculate {z{t)) and Xz(^) iii the regime 0(n~2/^) -C e -C 
0{n^), where 9 is expected to obey the Gaussian distribution 

P{9) = — e (54) 
with the normalization constant Zg. Defining the Fourier transform of Zu{t) as 

Zuit) = I (55) 
we write approximate expressions 

{z{t))c^ I ^5,(a;)e-*(e— ^), (56) 

* Some trajectories satisfy z{oo) > Zc, even when R < Re- Such a behavior can be described by (|43p. 
but not by ([50]) . Therefore, the argument below cannot be apphed to the calculation of (z(oo)), for 



example. Nevertheless, we expect that statistical properties of 6 are described by ([50]) when we restrict 
the trajectories z{oo) ~ 0. 
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and 

{4tr) ^J^J %Uu^)U^'y'-^-''' (e-(--')^) . (57) 
The Gaussian distribution fl54|) immediately leads us to 

{z{t)) I ^5.(.;)e^(*-<^))e-^, (58) 

and 

(zitr) ^J^J |^5„(..)5.(.')e'(-^-'^^*-<^»e-^^. (59) 
Using these, we can derive 




Here, let Tw = \/xe/n be the width of the distribution of 6. We expect that {z{t)) 
can be estimated as z^it — {9)) in the regime Tw <^ 1, where this regime is expressed as 
0(n-2/5) < e < 0(n°). Then, from we obtain 

X.(t) ~ Xe{dtz^{t))\ (61) 

from which we find that Xz(^) takes a maximum at t = t*, where (^^) ~ e^^/^ and 
X2(t=„) ~ — e~^/^. Since the fluctuation intensity is defined as the value of Xzi^) 
in the limit n oo with small e fixed, we conclude that ( = 1/2 and 7 = 5/2. The 
behavior of x{t) the regime 0(n^^/^) ^ e ^ 0(n~^/^), which is described by ( 160|) . 
seems complicated. We conjecture that there is no power-law behavior in this regime. 



6. Concluding remarks 

We have demonstrated that the edge deletion processes of random graphs exhibit the 
saddle-node bifurcation in the deterministic limit, as shown in (HTl) . The discontinuous 
transition of (h^t = 00)), from to Q{zc), is understood from the nature of the 
bifurcation of trajectories of z{t). (See ( |35l) for the transformation from z to h in 
the deterministic description.) We can also understand the divergent behavior of x{t) 
on the basis of critical fluctuations of exit time from the marginal saddle associated with 
the saddle-node bifurcation. 

The numerical analysis of the power-law divergences is quite difficult, although the 
increasing trends of xi'^) ^ind r are easily observed, as already shown in Fig. [2l For 
example, consider the power-law divergence x{^) — e"^/^ in the regime 0.01 < e < 0.1. 
In this case, we need to investigate the system with n ^ 10^. However, since 
our computational algorithm does not involve any tactical steps, we cannot perform 
numerical experiments of such a large system. 

Nevertheless, in Fig. |U we present the numerical result of random graphs on tt, = 2^^ 
vertices. The square symbols represent r and %(''") for several values of e. In order to 
complement the numerical data, we also display the results of numerical simulations of 
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Figure 4. (a) t as a function of e. (b) xi^) a-s a function of e. The square symbols 
represent tlie numerical results of the k-coie percolation dynamics with n — 2^^. 
The solid and dotted curves represent the numerical results of a Langevin equation 
with a small noise intensity that corresponds to the cases n = 2^^ and n ~ 2^^/3, 
respectively, in the fc-core percolation problem. The guide lines represent r = 8e^^'^ 
and X — O.le^^-^. 

a simple Langevin equation whose local form near the marginal saddle is equivalent to 
§0^. The equation is dtcp = -0((0- 1)^ + e) + V2r^, where (^(i)^(t')) = 5{t-t') and T 
is the noise intensity which is expected to be proportional to 1/n in the present problem. 
The solid curve corresponds to the case T = 2^^^/3, whose value is chosen such that 
the square symbols are on the solid curves. Then, by decreasing the noise intensity to 
T = 2"^°, which might correspond to the fc-core problem with n ~ 2^^/3, we obtain 
the dotted curve. As discussed theoretically, the power-law behavior of x — e~^/^ is 
observed in the regime 0(n-2/5) < e < 0(n°), and r ~ e is observed in the regime 
0(n-2/3) < e < 0(n°). 

With regard to finite size effects, we mention that the probability of finding 
trajectories satisfying /i(oo) 7^ is given by a universal function of n^/^(e— 2n~^/'^)l3 This 
implies that the system behavior in the regime e < 0{n~^^'^) is qualitatively different 
from that in the regime e > 0(ra~^/^). Theoretically, in order to describe the crossover 
around 0(n~^/^), we need to analyze fH3|) . not (l50l) . When we are interested in the 
relaxation behavior, we should focus on the regime e > 0{n~^^'^). 

It is worthwhile to note that xo exhibits the simpler behavior than x(t). We 
therefore conjecture that xe is a more fundamental quantity than x{t)- We also mention 
that the critical behavior of exit time from a marginal saddle is observed in a coupled 
oscillator model |T7j related to neuronal avalanches [18]. (See also Ref. [I9].) It is an 
interesting subject to find other examples belonging to the same universality class. 

Although the fc-core percolation is not directly related to jamming transitions, our 
results might provide a suggestion for future studies on jamming transition. As one 
example of such a study, we may theoretically consider the numerical result obtained 
for the jamming transition in the Frederic- Andersen model in a random graph |20j . 

ji This result was confirmed numerically by the direct simulations of the dynamics we consider. See 
also Ref. [16j as a mathematical argument. 
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because the fc-core percolation dynamics is regarded as an irreversible version of a kinetic 
constraint model. As another direction of study, one may analyze fluctuations of exit 
time in more general jamming systems. The important example is the application to 
the spherical p-spin glass model, for which the mode coupling theory is believed to be 
exact [21]. Since the transition described by this theory is interpreted as a variant of 
saddle-node bifurcation [22], we might discuss the divergent behavior of the so-called 
nonlinear susceptibility X4 on the basis of the exit time from the plateau regime. 

Finally, we consider the fc-core percolation in finite dimensional systems. In general, 
one may conjecture that a transition is smeared in a manner similar to bootstrap 
percolation problems [23]. (See Refs. [21] and [25] for attempts of studying the k- 
core percolation in finite dimensional systems.) From our viewpoint, as the first stage of 
a study on finite dimensional systems, we should identify the upper-critical dimension 
dc for a diffusively coupled model of a simple stochastic system undergoing a saddle- 
node bifurcation. With regard to this problem, we point out that z/* might be related 
to the exponent u that characterizes the divergence of the length scale as u = u^/dc 
[26] . Furthermore, in the next stage of studying finite dimensional systems, we should 
consider an equivalence or inequivalence between such a coupled model and the fc-core 
percolation dynamics in a finite dimensional lattice. The analysis developed in the 
present study may be useful in this consideration. 
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